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Abstract 

Background: Computational simulation using numerical analysis methods can help 
to assess the complex biomechanical and functional characteristics of the mitral 
valve (MV) apparatus. It is important to correctly determine physical contact 
interaction between the MV apparatus components during computational MV 
evaluation. We hypothesize that leaflet-to-chordae contact interaction plays an 
important role in computational MV evaluation, specifically in quantitating the 
degree of leaflet coaptation directly related to the severity of mitral regurgitation 
(MR). In this study, we have performed dynamic finite element simulations of MV 
function with and without leaflet-to-chordae contact interaction, and determined the 
effect of leaflet-to-chordae contact interaction on the computational MV evaluation. 

Methods: Computational virtual MV models were created using the MV geometric 
data in a patient with normal MV without MR and another with pathologic MV with 
MR obtained from 3D echocardiography. Computational MV simulation with full 
contact interaction was specified to incorporate entire physically available contact 
interactions between the leaflets and chordae tendineae. Computational MV 
simulation without leaflet-to-chordae contact interaction was specified by defining 
the anterior and posterior leaflets as the only contact inclusion. 

Results: Without leaflet-to-chordae contact interaction, the computational MV 
simulations demonstrated physically unrealistic contact interactions between the 
leaflets and chordae. With leaflet-to-chordae contact interaction, the anterior 
marginal chordae retained the proper contact with the posterior leaflet during the 
entire systole. The size of the non-contact region in the simulation with leaflet-to- 
chordae contact interaction was much larger than for the simulation with only 
leaflet-to-leaflet contact. 

Conclusions: We have successfully demonstrated the effect of leaflet-to-chordae 
contact interaction on determining leaflet coaptation in computational dynamic MV 
evaluation. We found that physically realistic contact interactions between the leaflets 
and chordae should be considered to accurately quantitate leaflet coaptation for MV 
simulation. Computational evaluation of MV function that allows precise quantitation 
of leaflet coaptation has great potential to better quantitate the severity of MR. 

Keywords: Mitral valve, Finite element, Contact interaction. Coaptation, 
Echocardiography 
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Background 

The mitral valve (MV) has the most complex structural characteristics of the four hu- 
man heart valves [1-3]. All components of the MV apparatus (mitral annulus, anterior 
and posterior leaflets, chordae tendineae, and papillary muscles) have very complicated 
interrelations that work together to provide the optimized valvular function. Abnormalities 
in any of these components may lead to valvular dysfunction: the most common pathologies 
include prolapse, regurgitation, annular dilation, and chordal rupture. In regurgitation, 
a well-coordinated action of the interrelated components of the MV apparatus plays 
an important role in preventing systolic backflow of blood into the left atrium [4] . 

Advances in ultrasound imaging techniques have allowed clinicians to accurately 
assess MV physiology, and have aided in determining the degree of valvular ab- 
normalities and the necessity for appropriate surgical intervention. Conventional 
two-dimensional (2D) echocardiography provides excellent insight into the complex car- 
diac anatomy in real time [5]. Three-dimensional (3D) echocardiography provides further 
detailed morphology of the MV apparatus better contributing to our understanding of the 
anatomical and physiologic characteristics of MV function [6]. However, biomechanical 
information (e.g., high stress concentration, abnormal strain distribution, malformed leaf- 
let bending deformation, etc.) over the MV structure cannot be obtained from echocardi- 
ography. It is well known that mechanical stress is one of the primary triggering factors in 
deterioration and failure of heart valve tissue [7-10]. Computational simulation using nu- 
merical analysis methods can help to assess the complex biomechanical and functional 
characteristics of the MV apparatus [8,9,11-16]. Finite element (FE) analysis has been 
employed to determine deformation patterns and stress distributions across heart valve 
tissue under the physiologic conditions [11]. 

Using FE analysis and computational simulation, it is important to correctly model 
the contact interaction between the MV apparatus components. We have recently dem- 
onstrated that failure of proper leaflet coaptation leads to reversal of blood flow from 
the left ventricle into the left atrium, i.e., mitral regurgitation (MR) [17,18]. Contact 
interaction between the two MV leaflets is generally modeled by surface-to-surface 
contact such that two master and slave pairs are enforced by the penalty constraint to 
characterize appropriate contact interaction between the soft and wet leaflet surfaces 
[8,15,16,19-21]. There is an additional type of primary physical contact interaction between 
the MV apparatus components: leaflet-to-chordae contact interaction. A proper computa- 
tional modeling strategy for this contact interaction between the leaflets and chordae tendi- 
neae is required to more accurately perform simulations of MV function. Previous 
computational MV studies focused primarily on the contact interaction between the MV 
leaflets as an interface pair [8,15,16,19-21]. However, the effect of leaflet-to-chordae contact 
interaction on computational MV evaluation has not been ftiUy investigated. 

We hypothesize that leaflet-to-chordae contact interaction plays an important role in 
computational MV evaluation, specifically in quantitating the degree of leaflet coapta- 
tion directly related to the severity of MR. 

In this study, we have acquired two different MV types from patients using 3D echo- 
cardiographic data, modeled virtual MVs, performed dynamic FE simulations of MV 
function with and without incorporation of leaflet-to-chordae contact interaction, and 
determined the effect of leaflet-to-chordae contact interaction on the computational 
MV evaluation. 
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Methods 

Virtual MV modeling 

The Committee for the Protection of Human Subjects at The University of Texas 
Health Science Center at Houston approved our translational research protocol to 
utilize 3D transesophageal echocardiography (TEE) for computational modeling of pa- 
tient MVs. A Philips iE33 ultrasound system (Philips Medical Systems, Bothell, WA) 
with a 3D TEE transducer was utilized for imaging of the patient MV apparatus. Fol- 
lowing our standard clinical protocols, 3D TEE was performed to acquire the MV geo- 
metric data in two patients. One patient had a normal MV without MR, and another 
had a pathologic MV with MR. MV modeling was conducted by our virtual MV model- 
ing protocol [17,18]. Briefly, the geometric information of the MV apparatus including 
annulus, leaflets and papillary muscles at end diastole (open) was identified from 3D 
TEE data, segmented, and traced in the cylindrical coordinate system using a custom- 
designed semi-automated image processing algorithm. The geometric data was trans- 
formed into the Cartesian coordinate system, followed by virtual MV modeling. The 
mitral leaflets and annulus at end diastole were created utilizing the non-uniform ra- 
tional B-spline (NURBS) surface modeling technique, and meshed with 3D triangular 
shell elements. The chordae tendineae were modeled by adding 3D line (truss) ele- 
ments between the papillary muscles and the leaflet edge. 

Material modeling and boundary conditions 

Detailed protocols of material modeling and dynamic FE simulation of MV function 
are described in a previous study [18]. The MV leaflet tissue was modeled as an aniso- 
tropic hyperelastic material using a general Fung-type elastic constitutive model. The 
stress-strain relationship of the leaflet tissue was defined along the circumferential and 
radial directions. Material parameters were determined by fitting the constitutive model 
to the biaxial mechanical test data of the anterior and posterior leaflet tissue from a 
previous study [22]. These experimentally determined material characteristics data of 
the anterior and posterior leaflet tissue were implemented into the Fung-type elastic 
material model in ABAQUS (SIMULIA, Providence, RI). Leaflet thickness was set to be 
0.69 mm and 0.51 mm for the anterior and posterior leaflets, respectively [23]. 

The chordae tendineae were modeled as nonlinear hyperelastic materials utilizing 
the first and second order Ogden models. Experimentally determined material charac- 
teristics of the posterior marginal chordae, anterior marginal chordae, and strut 
chordae from a previous study [19] were incorporated into our MV modeling, re- 
spectively. Cross-sectional areas were set 0.27 mm^ for the posterior marginal chordae, 
0.29 mm^ for the anterior marginal chordae, and 0.61 mm^ for the strut chordae 
[19]. Density and Poissons ratio of the tissue were set at 1,100 kg/m^ and 0.48, re- 
spectively [14,21]. 

Mitral annulus and papillary muscle configurations at end diastole (open) and peak 
systole (closed) were pre-determined in the 3D TEE data. Combined with time-varying 
nonlinear nodal displacements between end diastole and peak systole, dynamic bound- 
ary motion was defined for both the annulus and the papillary muscles [18]. Free rota- 
tional movement was allowed at the annulus and papillary muscle tips. Time-varying 
physiologic transvalvular pressure loading was applied onto the leaflets for dynamic 
simulation of MV function across the cardiac cycle. 
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Contact modeling 

The contact interaction in computational MV simulations primarily consists of leaflet- 
to-leaflet and leaflet-to-chordae contacts. In ABAQUS, the general contact algorithm 
allows to define edge-into-edge contact penetrations [24], The general contact algo- 
rithm creates contact forces to resist edge-into-edge penetration (i.e., interaction be- 
tween the shell perimeter edges and/or the chordae truss edges). The penalty method 
is used to enforce contact constraints in the general contact algorithm. The penalty 
contact algorithm searches for slave node penetrations and edge-to- edge penetrations, 
and equal and opposite forces are applied to the slave nodes to oppose the penetration 
[24]. This process involves the addition of penalty terms in FE analysis. 

Figure 1 demonstrates a contact pair with frictional contact, which has a slave node 
(point S) and a master line (connecting node 1 and 2) in the general contact algorithm. 
So and S refer to the slave node before and after loading, respectively [25]. The penalty 
method utilizes the penalty term which contains the predetermined stiffness matrix 
(penalty stiffness) of the contact surface. The stiffness matrix of the contact surface is 
incorporated into the stiffness matrix of the contacting body: 



(1) 



where is the stiffness matrix of contacting body, is the stiffness matrix of contact 
surface, u is the displacement, and F is the applied force [25]. Any violation of the con- 
tact condition is punished by increase of the total virtual work by the penalty term. 
The first variation of the potential energy of the contact element is 



(2) 



where kn is the penalty term along the normal direction, kt is the penalty term along 
the tangential direction, and gn and gt are the penetration and sliding along the normal 




fn 



Figure 1 A schematic of a contact pair (slave and master surfaces) and variables used in the 
penalty method. 
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and tangential directions, respectively [25]. Hence, the equation can be written as 
follows: 

fn^gn +ft^gt = ^ngnSgn + ^%^{gt)t^d^ngnSgt (3) 

fn = kngn (4) 

ft = -^%^{gt)t^d{kngn) (5) 

C -1 forg,<() 

sgn(g,) = 1 0 forg,=0 (6) 
I 1 forg,>0 

with a coefficient of kinetic friction. 

A schematic of the MV contact configuration under the normal physiologic conditions 
is displayed in Figure 2A. Physiologically appropriate leaflet-to-leaflet and leaflet-to- 
chordae contacts during systole are demonstrated. If FE analysis employs leaflet-to-leaflet 
contact alone (i.e., without incorporation of leaflet-to-chordae contact), penetration of 
the chordae tendineae into the leaflet structure is expected during the entire systolic 
phase (Figure 2B). 

Computational MV simulation with full contact interaction was specified such that 
the contact domain of the general contact algorithm was defined to incorporate entire 
available contact interactions between the leaflets and chordae tendineae. Computa- 
tional MV simulation without leaflet-to-chordae contact interaction was specified by 
defining the anterior and posterior leaflets as the only contact inclusion in the general 
contact domain. A frictional coefficient value of 0.05 was utilized to describe the fric- 
tional behavior between the contact tissue [20]. 



Results 

MV leaflet deformation 

Figure 3 exhibits leaflet deformation at peak systole with and without incorporation of 
leaflet-to-chordae contact interaction. The difference in leaflet deformation between 




(A) With Leaflet-to-chordae Contact (B) Without Leaflet-to-chordae Contact 

Figure 2 Contact configurations between the MV apparatus components with and without 
consideration of leaflet-to-chordae contact interaction. 

V ) 
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the two contact modeling methods was best demonstrated in the vicinity of the com- 
missural regions. The posterior leaflets of both normal and pathologic MVs were longer 
than the anterior leaflet in the commissural regions. With leaflet-to-chordae contact 
interaction incorporated, the anterior marginal chordae retained the proper contact with 
the posterior leaflet during the entire systole (Figure 3A). On the other hand, the anterior 
marginal chordae demonstrated unrealistic penetration into the posterior leaflet during 
systole without leaflet- to-chordae contact interaction incorporated (Figure 3B). In particu- 
lar, the pathologic MV with MR demonstrated a greater degree of chordal penetration in 
the anterolateral commissural region compared to the normal MV. 

Alteration in leaflet deformation due to the incorporation of leaflet-to-chordae 
contact interaction revealed a difference in leaflet coaptation of the pathologic MV 
(Figure 4). With leaflet-to-chordae contact interaction, the cross-sectional view of the 




(A) With Leaflet-to-chordae Contact (B) Without Leaflet-to-chordae Contact 



Figure 4 Leaflet coaptation of the pathologic MV involving MR at peak systole with and without 
incorporation of leaflet-to-chordae contact interaction. 

V J 
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leaflet deformation near the anterolateral commissural region at peak systole demon- 
strated incomplete coaptation between the anterior and posterior leaflets. This lack of 
leaflet coaptation disappeared when only leaflet-to-leaflet contact interaction was included 
(i.e., without leaflet-to-chordae contact). The maximum gap distance was 3.2 mm and 
1.4 mm with and without leaflet-to-chordae contact interaction, respectively. 



Leaflet contact distribution 

The extent of leaflet coaptation, which directly correlates with the degree of MR, was 
determined by measuring the distribution of contact pressure between the anterior and 
posterior leaflets at peak systole (Figure 5). A threshold value (50 kPa) of contact pres- 
sure upon the anterior leaflet was employed to best demonstrate fuU contact leaflet co- 
aptation. The maximum leaflet contact pressure decreased from 77 kPa to 62 kPa in 
the normal MV when leaflet-to-chordae contact interaction was incorporated. The dif- 
ference of maximum leaflet contact pressure in the pathologic MV was much greater 
(87 kPa vs. 52 kPa) between the simulations with and without leaflet-to-chordae con- 
tact interaction. 

The normal MV showed sufficient leaflet coaptation and little difference in overall 
leaflet contact pressure distribution between the simulations with and without incorp- 
oration of leaflet-to-chordae contact interaction. In the pathologic MV with MR, there 
was a large non-contact region near the anterolateral commissure, indicating a consid- 
erable occurrence of MR. The length of leaflet marginal edge with no coaptation in the 
anterolateral region was much larger in the simulation with leaflet-to-chordae contact 




Non-contact Region 
Contact Region 
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(A) With Leaflet-to-chordae Contact (B) Without Leaflet-to-chordae Contact 

Figure 5 Distribution of contact pressure between the anterior and posterior MV leaflets at peak 
systole with and without incorporation of leaflet-to-chordae contact interaction. 
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interaction (6.9 mm) than for the simulation with leaflet-to-leaflet contact alone 
(2.6 mm). This data corresponds to the markedly enlarged size of regurgitant orifice 
area (10.2 mm^ vs. 1.4 mm^) when leaflet- to-chordae contact interaction was imple- 
mented. With incorporation of leaflet-to-chordae contact interaction incorporated, the 
size of total contact region decreased by 3.8% and 4.2% in the normal and the patho- 
logic MVs, respectively. This indicates that leaflet-to-chordae contact interaction affects 
both qualitative and quantitative determination of leaflet coaptation when evaluating 
MV function. 

Leaflet stress distribution 

The maximum principal stress distributions across the mitral annulus and leaflets at 
peak systole are demonstrated in Figure 6. The normal and pathologic MVs showed a 
difference in stress distribution pattern. The maximum stress values in the normal MV 
occurred near the anterior saddle horn region spreading along the circumferential dir- 
ection. In the pathologic MV, a much wider range of large stress distribution was found 
across both leaflets spreading along the radial direction. This directional pattern of 
stress distribution was clearly observed around the regurgitant area of the anterior leaf- 
let close to the anterolateral commissural region. The maximum stress values in the 
pathologic MV was twice larger than for the normal MV. 

MV simulations with and without consideration of leaflet-to-chordae contact inter- 
action demonstrated little difference in overall stress distribution in both normal and 
pathologic MVs. The maximum stress value (0.32 MPa) at peak systole was comparable 
in the normal MV simulations with and without leaflet-to-chordae contact interaction. 




Pathologic MV 

(A) With Leaflet-to-chordae Contact (B) Without Leaflet-to-chordae Contact 

Figure 6 Leaflet stress distribution with and without Incorporation of leaflet-to-chordae 
contact Interaction. 
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However, in the pathologic MV, the maximum stress value decreased from 0.7 MPa to 
0.66 MPa when leaflet-to-chordae contact interaction was incorporated. A wider range 
of large stress concentration was found in the regurgitant areas beside both leaflets. 

Discussion 

Accurate determination of the extent and severity of MR is important as MR can cause 
serious pathophysiologic problems leading to left ventricular dysfunction and heart fail- 
ure. The most popular modality for MR evaluation at present is color Doppler ultra- 
sound imaging that provides semi-quantitative information of MR occurrence [26]. 
Computational simulations of MV function with the use of patient-specific MV geomet- 
ric data obtained from clinical imaging modalities can provide additional biomechanical 
and functional information [17,18,27,28]. Biomechanical information such as leaflet con- 
tact pressure distribution and the extent of leaflet coaptation can directly provide quanti- 
tative information pertaining to the degree and mechanism of MR [17,18]. Therefore, it is 
important to correctly determine physical contact interaction between the MV apparatus 
components during computational MV evaluation. 

In the present study, we investigated the effect of leaflet-to-chordae contact inter- 
action on patient-specific computational MV evaluation. To the best of our knowledge, 
the effect of leaflet-to-chordae contact interaction on computational MV evaluation has 
never been fully investigated. 

The greatest difference between the computational MV evaluations with and without 
leaflet-to-chordae contact interaction was found in the deformation of the leaflets and 
chordae as well as in the distribution of leaflet coaptation. Without incorporation of 
leaflet-to-chordae contact interaction (i.e., consideration of leaflet-to-leaflet contact only), 
the computational MV simulations demonstrated physically unrealistic contact interactions 
between the leaflets and chordae (Figure 3). The chordal structures penetrated into the leaf- 
lets as soon as contact between the two structures occurred during the systolic phase. It is 
clear that this non-physical geometric alteration can overestimate leaflet-to-leaflet contact 
and lead to inaccurate determination of leaflet coaptation (Figure 5). This increase in leaflet 
coaptation in computational MV evaluation can directly translate into underestima- 
tion of the severity of MR [17,18]. Several studies recently reported computational evalu- 
ation of MV function accompanied with patient-specific MV geometric data obtained 
from clinical imaging modalities [17,18,27,28]. The patient MVs in these studies demon- 
strated a high degree of asymmetric structural characteristics in the mitral annulus, 
leaflets and chordae distribution. The present study demonstrates that a disregard of 
leaflet-to-chordae contact interaction can affect computational evaluation of MV fiinction, 
particularly in pathologic MVs involving considerable MR. 

Stress distribution was evaluated for the normal and pathologic MVs to compare the 
effect of leaflet-to-chordae contact interaction on alteration of biomechanical charac- 
teristics of the MV apparatus structure during valve function (Figure 6). The pathologic 
MV revealed an apparent stress distribution pattern along the radial direction across 
both the leaflets while the normal MV showed a relatively uniform stress distribution. We 
speculate that this particular stress distribution pattern in the pathologic MV was induced 
by the dilated mitral annulus and the increased tensile stresses in the adjacent tightened 
chordae in the region where MR occurred. Incorporation of leaflet-to-chordae contact inter- 
action revealed a relatively small difference in the overall pattern and magnitude of stress 
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distribution compared to those without the contact interaction. This may be the reason 
why previous computational MV evaluation studies did not pay attention to leaflet-to- 
chordae contact interaction, and described little details on contact interaction modeling 
between the interrelated components of the MV apparatus [19,21,29]. We found that the 
difference in stress distribution may be negligible during valve function in normal MVs 
regardless of clear incorporation of leaflet- to -chordae contact interaction. However, there 
is great potential to generate a considerable error in determining the maximum stress 
concentration and stress distribution in pathologic MVs when leaflet- to- chordae contact 
interaction is not correctly incorporated in the computational simulation. 

Although the present study has demonstrated the effect of contact interactions be- 
tween the MV apparatus components on computational MV evaluation, the simulation 
results should be carefully interpreted prior to defining MV pathophysiology. Although 
we demonstrated computational MV evaluation using two patient data (one normal 
and one pathologic MVs), the physical influence of leaflet- to- chordae contact inter- 
action should be effective in any types of MV simulation studies. We utilized previously 
published material characteristic data, employed geometric and physiologic parameters 
such as leaflet thickness and tissue density, and assumed homogeneous thickness and 
material behavior of the MV tissue. Although different types of tissue material behavior 
may affect computational MV simulation outcomes, it would not change the primary 
findings of this study in terms of the physical effect of leaflet- to- chordae contact inter- 
action. Our computational simulation data was obtained purely by structural ana- 
lysis, i.e. dynamic FE simulation incorporated with physiologic damping conditions. 
The latter implies that blood fluid effects on normal and shear pressure distribution over 
the MV leaflets and chordae were not fully incorporated in the computational simula- 
tions. Notwithstanding these limitations, the effect of leaflet-to-chordae contact inter- 
action on computational MV evaluation presented in this study should be valid in any 
types of MV simulation studies. 

Conclusions 

We have successfully demonstrated the effect of leaflet-to-chordae contact interaction 
on determining leaflet coaptation in computational dynamic MV evaluation. Two types 
of patient-specific MV geometric data were analyzed to assess the difference of this model- 
ing effect in normal and pathologic MVs. We found that physically realistic contact interac- 
tions between the leaflets and chordae should be considered to accurately quantitate leaflet 
coaptation for MV simulation. Computational evaluation of MV function that allows precise 
quantitation of leaflet coaptation has great potential to better quantitate the severity of MR. 
Improved MR quantitation using this computational simulation combined with patient- 
specific 3D echocardiography can help us to better diagnose and treat MV pathology. 
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